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In a suitably chosen back-to-back kinematics, four-jet production in hadronic collisions is known 
to be dominated by contributions from two independent partonic scattering processes, thus giving 
experimental access to the structure of generalized two-parton distributions (2GPDS). Here, we 
show that a combined measurement of the double hard four-jet cross section in proton-proton and 
proton-nucleus collisions will allow one to disentangle different sources of two-parton correlations in 
the proton, that cannot be disentangled with 4-jet measurements in proton-proton collisions alone. 
To this end, we analyze in detail the structure of 2GPDS in the nucleus (A), we calculate in the 
independent nucleon approximation all contributions to the double hard four-jet cross section in 
pA, and we determine corrections arising from the nuclear dependence of single parton distribution 
functions. We then outline an experimental strategy for determining the longitudinal two-parton 
correlations in the proton. 
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I. INTRODUCTION 



With increasing center of mass energy at hadron colliders, multi-particle final states at high transverse momentum 
receive an increasingly important contribution from multi-parton interactions. The prototype of such processes are 
double parton interactions in which two partons from each hadron enter into collision at two distinct hard vertexes. 
The calculation of such double parton interactions involves double Generalized Parton Distributions (2GPDS) that - in 
principle - contain information about the spatial and momentum correlations of two partons in the incoming hadronic 
wave function. In practice, experimental constraints on these correlations are scarce. Since nuclear projectiles offer 
significantly different two-parton correlations, the question arises to what extent a combined analysis program of 
proton-proton and proton- nucleus collisions at the LHC could help to constrain 2GPDS. In this paper, we classify the 
different contributions to nuclear 2 GPDs and how they contribute to the double parton scattering cross section into 
4 jets that is experimentally accessible in the upcoming p-A run at the LHC. 

We focus on the physics in a suitably chosen back-to-back kinematics: the dijet momentum imbalances <5 2 3 , <5| 4 are 

^13 = (&it + fat) 2 ^ kit ~ ^It ~ *9i> ^24 = {fat + kit) 2 -C fc 2 t ~ feft ~ Qis where the momenta kit, fat, fat, fat are the 
transverse momenta of individual jets, the scales Qi,Q% are the resolution scales of two partons. In this kinematics the 
production of four jets in the collision of two partons (so-called 2 — > 4 process) is suppressed in the leading logarithmic 
approximation compared to the hard processes that involve the collision of four partons [IH5] . For this back-to-back 
kinematics, it is useful to write 2GPDS as the sum of two contributions, 

2 G p ( Xl ,x 2 , Ql Ql A) = G«° Me (x 1 ,x 2 , Ql, Ql A) + Gf^ c (x 1 ,x 2 ,Q\ 1 Q 2 , A) . (1) 

Here, 2GPDS are written as functions of the momentum fractions x\,x 2 and the resolution scales Q\,Q 2 of the two 
partons. In the following, we will not write explicitly the dependence of 2GPDS and GPDs on resolution scales. 
The transverse momentum parameter A denotes the difference between the transverse momenta exchanged by one 
parton in the amplitude and complex conjugate amplitude and it is conjugate to the relative spatial transverse 
distance between the two partons. In general, 2 GPDs contain a non-perturbative two-parton contribution in which 
QCD evolution amounts to an independent evolution of both partons with standard one-parton evolution equations 
(Gp° ubl °). In addition, there is a contribution in which both partons result perturbatively as the two daughters of a 
single parent parton in the QCD evolution (Gp lngl0 ). This second term couples the evolution of 2 GPDs to the evolution 
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of the standard one-parton distribution functions. For the production of two pairs of jets in independent back-to-back 
kinematics, both contributions are known to contribute with parametrically equal weights [3j. 

The double hard four jet cross section for the collision of hadrons A and B can then be written in terms of 2 GPDs 

as 

d(J i 3 B et f d 2 A d&xix'^xx) da 2 (x' 2 ,x 2 ) n , , , n . x , , . 

' " ■ 2 G A (x 1 ,x 2 , A) 2 G B (x 1 ,x 2 , A) . (2) 



dt x dt 2 J (2?r) 2 dt x dt 2 

We denote by &i the partonic 2 — > 2 scattering cross sections. The measurable four-jet cross section in AB collisions, 
is a function of the four jets' cm. transverse energy and rapidity, that are connected to the variables Xi, x[ 
and the virtualities in the standard way. For the case that <7i and o 2 denote indistinguishable scattering processes, 
equation ^ must be multiplied by a symmetry factor 1/2 that we omit for briefness. In proton-proton collisions, it 
is customary to parametrize the cross section o™ et as the product of two two-jet cross sections o" 2 j et , 

daf oet _ 1 daf jet daf jet 
di\di 2 S dtx dt 2 

Here the quantity S (sometimes referred to as <7 e ff) characterizes the effective transverse area of the four parton 
interaction and the effect of longitudinal correlations between the partons in the colliding hadrons. In general, S can 
depend on the momentum fractions Xi and virtualities of the incoming partons. Data from the Tevatron indicate that S 
is typically of the order of 15 mb [IHZ]. This is a factor ~ 2 smaller than expectations based on uncorrelated two-parton 
distributions, see e.g. [5]. It is a clear indication that non-trivial parton correlations are experimentally accessible in 
double hard parton interactions. To date, it remains an open question of whether these two-parton correlations in 
the proton are predominantly transverse (as realized e.g. in models that picture the proton as composed of several 
hot spots) or predominantly longitudinal (as resulting e.g. from perturbative 1 — > 2 splittings). These and other 
questions about the Tevatron data have motivated much work recently [M2"5] . Older relevant work includes [STJ . 

The case of double hard four jet production in proton-nucleus collisions was discussed first in |28j where it was 
pointed out that such measurements would be sensitive to the longitudinal correlations of the partons. More recently, 
multiple parton interactions were considered in the production of two leading pions in deuteron-gold collisions at 
RHIC 29J, and for the case of proton-deuteron collisions [3D]. In the present paper, we aim to extend this discussion 
in the light of the recent pQCD studies [THS] mentioned above. For the case that A is a nucleus, the 2GPD may be 
written as the sum of three distinct contributions, depending on whether both partons belong to the same nucleon 
(IN) or to different nucleons (2N). Distinguishing for the first case again single from double contributions, as in eq. ([!]), 
one obtains 

2 G A ( Xl ,x 2 , A) = G A nslc ' 1N (x!,x 2 , A) + G d A mhie > m (xx,x 2 ,A) + G 2 f(xi,x 2 , A) . (4) 

As will become clear in the following, this classification relies on viewing the nucleus as a superposition of nucleons 
to which partons can be asigned uniquely. Inserting the expressions for 2GPDS inside a nucleus Q and inside a 
nucleon ([!]) into the double hard four jet cross section |2]), the following different contribution emerge for the case of 
proton-nucleus scattering: 

T ^double (gj ^.double, IN 



II. 


/^double 
U p 


® G A 


t 


le,lN 


III. 


/^double 
Up 


®G\ N 


direct 


IV. 


G p inglc ( 


E> Gf g 


le,lN 


V. 


G si„glo ( 


y-ydouble,lN 


VI. 


G p inglc ( 


*>G\ N 


direct 



The terms in which two partons are taken from the same nucleon in a nucleus contribute equally to proton-nucleon 
collisions. More specifically these are the 4 — > 4 contribution (case I) and the 3 — > 4 contributions (case II and V). 
They will be discussed in section [TTJ It is known that the single-single contribution (case IV) does not make a dominant 
contribution to the back-to-back kinematics; rather, it can be regarded as a one-loop correction to the 2 — > 4 cross 
section [T]. 




FIG. 1: Schematic representation of the contributions G^ nglc ' 1N + G A ouhlc ' 1N to the nucleus 2GPD that enter the cross section 
for double hard 4-jet production (cases I, II and IV, V discussed in the text). Two patrons of momentum fractions Xi, X2 
are drawn from the same single nucleon in the amplitude (left hand side of the diagram) and complex conjugate amplitude. 
Depending on whether the two nucleons arise perturbatively in a 1 — > 2 splitting from a single parton in this nucleon, or 
whether they are of non-perturbative origin, we shall refer to them as G^ ngle ' 1N and G A ouhlc ' 1N , respectively. 



For a nuclear projectile, two additional contributions arise. These involve two partons from different nucleons in the 
nucleus that interact with Gp° uUe (case III) or G p mslc (case VI) of the proton, respectively. These will be discussed 



in section III We have labeled both these contributions with the subscript 'direct' to indicate that the partons with 
momenta x\ and x 2 are assigned to the same nucleons in amplitude and complex conjugate amplitude. 

Viewing the nucleus as a bound state of many nucleons (without considering modifications to their internal structure) 
is an approximation. If we relax the working hypothesis that partons can be assigned uniquely to nucleons in a nucleus, 
then additional contributions are possible. First, it is possible that two nucleons of the nuclear wave function are 
involved in both amplitude and complex conjugate amplitude, but that the two partons are interchanged across the 
cut. We label these contributions with the subscript 'interference'. Second, it is conceivable that the two partons 
taken from the nucleus belong to one nucleon in the amplitude but to two different nucleons in the complex conjugate 
amplitude. 
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IX. 1N-2N interference 

In general, such interference terms indicate the break-down of a probabilistic picture of double-parton interactions in 
pA collisions in terms of an independent superposition of nucleons. By their very nature, they characterize partonic 
cross-talk between different nucleons in a nucleus and may thus help to elucidate partonic nuclear structure. We note 
that deviations from the simple picture of a nucleus as a superposition of nucleons are known already on the level of 
single parton distributions, e.g. as EMC and nuclear shadowing effects [31]. In the following sections, we discuss the 
different contributions to the double-hard four jet cross section following the classification listed above. 



II. SINGLE NUCLEON SCATTERING TERMS (CASES I, II AND V) 

For the contributions I. -VI., the momenta of the nucleons are the same in amplitude and complex conjugate 
amplitude. Therefore, if both partons belong to the same nucleon (cases I, II and V), one can integrate over the 
momenta of all other nucleons and write the corresponding part of the nuclear 2GPD as 

/-<1N 1 \\ / 1 /Visingle/^l x 2 A /-.double/ x l x 2 a\\ N I j2 /r\ 

G a (aJi,x 2) A)= / -J [G N B — ,— ,A) + G N ( — ,—,A))p A (a,p t ) — dp t . 5) 
J or V a a a a ) a 

Here, the nuclear 2GPD is the sum of the terms Q s ™9 le -- 1N anc [ G A ouhlc ' 1N , introduced before, and expressed as 
an integral over the corresponding contributions of individual nucleons. The quantity p A (a,Pt) denotes the light- 
cone nucleon density of the nucleus normalized as J p A (a,pt)da/a = A. The factor l/a for each of the partons 
reflects the fact that the number of partons between x± and X2 should not change under Lorentz boosts. These 
extra factors l/a are absorbed in the flux such that one recovers for the structure function the standard expression 
F 2A (x,Q 2 ) = J p N A {a, Pt )^d 2 p t F 2N ^Q 2 ). 




FIG. 2: Schematic representation of the contribution G \ to the nucleus 2GPD that enters the cross section for the double hard 
4-jet production in the terms III and VI discussed in the text. Here, the two patrons are drawn from two different nucleons 
that carry longitudinal nucleon momentum fractions ot\ and 02- 



Equation ^ is written for an ensemble of A moving nucleons satisfying the momentum sum rule 
J ap^(a,pt)da/ad 2 p t = A since a « = A. This raises the question of how well one can approximate G 1 ^ 1 in 
terms of the distributions G s ™ 9le (xi, x 2 , A) + Gj^ uhlc (xi, x 2 , A) written in the nucleus rest frame. To address this 
question, we replace a = 1 + (a — 1) in the arguments of (1/a 2 ) G 1N — , A), and we expand in powers of (a — 1). 
Using momentum sum rule and baryon number sum rule, we find J —(a — l)p^(a) — A — A — 0, see Ref. [521 133] . 
Therefore, corrections due to Fermi motion (i.e., corrections due to the a-dependence of the integrand of G 1 ^) arise 
only to second order in (a — 1). The longitudinal momentum distribution of nucleons in a nucleus peaks at a = 1 
with small dispersion, and therefore 

G 1 /(x 1 ,x 2 ,A)=AG N (x 1 ,x 2 ,A)(l + On (a - \) 2 p N A { a , Pt ) d ^d 2 p^ . (6) 

Here, the correction term involves first and second derivatives of the nucleon 2GPD with respect to x\ and x 2 . A 
precise numerical estimate will have to constrain this term numerically. Parametrically, the correction is small. The 
dominant linear dependence of G 1 A N (xi, X2, A) on nucleon number A translates directly into a linear dependence of 
the corresponding contribution to the double hard four-jet cross section 

~ A da ^ e l = A d<T ^ et da ™ et . (7) 
dti dt 2 dti dti & dti dt 2 

^pAAN 

We have introduced an m sign in this relation to indicate that the identification of .p't- with A times the cor- 

dti dt 2 

responding cross section in pp relies on neglecting the nuclear modification of parton distribution functions. Our 
discussion up to section VI will rely on this approximation. This is justified since we consider larger Q 2 processes 



at moderate x, where nuclear modifications are expected to be small. Within this approximation, the contributions 
to the double hard four jet cross section in pA discussed here are exactly the contributions that one obtains from 
superimposing four jet cross sections from A independent nucleon-nucleon collisions; the effective transverse area S 
in |7]) is therefore the quantity measured in pp collisions. In section VII we go beyond this approximation and we 



discuss how the nuclear dependence of parton distribution functions can be taken into account. 



III. DOUBLE NUCLEON SCATTERING TERMS (CASES III AND VI) 



Figure |2| shows the 2GPD contribution in which the two partons belong to two different nucleons in both 
amplitude and complex conjugate amplitude. This term enters the contributions III and VI of the double hard four 
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jet cross section. In terms of the nuclear light cone wave function ip A of the A-nucleon system, it takes the form 



G A "( Xl ,x 2 ,A) = A(A-l) 



i=A 

-n 



daid 2 p ti 



a x a 2 



^a 4 -A 5 (2) {^pii ip A (ai,a 2 ,p n ,p t2 ,-.) (8) 



xip A (a-L,a 2 ,Pti + A,p t2 - A, ....)G N {xi/a u |A|)GAr(a;2/a2, |A|). 

Here, the transverse momentum transfer A is exchanged between the two active nucleons, while their light cone 
fractions are conserved. The partonic momentum fractions Xi drawn from the two active nucleons are then determined 
via the generalized parton distributions Gjsi{xi/ai, A) of the nucleons. Since the wave function ipA is normalized to 
unity, the prefactor A (A — 1) results from combinatorics. The factor 1/ (ot\ a 2 ) has the same origin as the factor I /a 2 
in eq. 

Expanding in equation ^ the arguments of Gat in powers of [oti — 1), the leading term can be written in a factorized 
form involving the two-nucleon form factor i^ oubl ° 



G A N ( Xl ,x 2 ,A) = A(A-l)G N ( Xl ,\A\)G N (x 2 ,\A\) i^ ouble (A, -A) , 



double / 



redouble 



(A, -A) 



n daid 2 pu 



\Y^ a *~ A J S(2) ^5Z Pti j ^A( a l: a 2,Ptl,Pt2,----) 

xip A (ax,a 2 ,p t i + A,p t2 - A, ....) . 



(9) 



(10) 



As for subleading correction, we note that one can exploit the symmetry of the integrand of ([£]) under a, — >(2 — a.i) 
in the nonrelativistic limit to see that corrections to M are of the order (1 — a) 2 , similar to the case of eq. (pi). Since 
the momentum fraction on of all nucleons in the nucleus are close to unity, we can approximate them in the nucleus 
rest frame in the non-relativistic limit, at = 1 + psi/m^- The two-nucleon form factor reads then 



F 



NR double 



A 



(A, —A) = J fj] d 3 P )j r A (PuP2, ...)^ A {px + A, P2 - A,p 3 , ...)S^ (j^pi 



(11) 



Equations (10) and (11 ) express the two-nucleon form factor for an arbitrary nucleus wave function, and therefore can 
account for arbitrary nucleon correlations. Also, when combined with eq. ([£]), these expressions can account for the 
finite size of nucleons as well (which is characterized by the single nucleon GPDs). We now discuss how more explicit 
expressions, suitable for direct numerical evaluation, can be obtained if assumptions about nucleon correlations and 
the finite size of nucleons are made. 

First, we turn to the independent nucleon approximation, when the nuclear wave function is written as a product 
of single nucleon wave functions. This neglects all internucleon correlations, including constraints from recoil that 
arise from the kinematic 5-function in ( 11 ). (These latter corrections are proportional to 1/A. A parametrically more 
important source of corrections to this picture of double hard 4- je t pr oduction arises from short-range NN interaction 
that are suppressed by a factor oc 1/A 1 / 3 .) One can express (11) in terms of products of Fourier transforms of 
single nucleon wave functions ipN^i)- Using a single nucleon density p A (r) = A ip%(r) ipN (r), that is normalized to 
J p A (r)d?r = A, the two-nucleon form factor is the product of single nucleon form factors 



redouble 
*A 



(A, 



d 3 r -^pA{r) exp 



iA ■ r 



F A (Af 



(12) 



Since A is a two-dimensional vector in the transverse plane, the single nucleon form factor can be written in terms of 
the nuclear thickness function T(b) — dzp A (z, b) as 



Fa(A) 



d 2 bT(b)exp{iAb). 



(13) 



The well-known approximate relation between the form factor and the nucleus radius Ra, 

1 



TT'doublc 
b A 



(A, -A) a exp 



-A 2 R 2 



(14) 



can then be obtained by expanding (|13| for small A, F A (A) ~ 1 — ^A 2 R 2 A and reexponentiating this expression. 
However, the Gaussian approximation (|14h somewhat underestimates the drop of F A (A 2 ) with A 2 for A 2 R\/6 > 1. 
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So, while ( 14 1 is well-suited for parametric arguments, it is preferable to base numerical estimates on evaluating ( 13 ) 
without further approximation. 

Second, we discuss now approximations that amount to neglecting the nucleon size in comparison to the nucleus 
size. According to (fl4| ), the double nucleon scattering contribution (J9j) to the 2GPD has its main support for small 
values of A 2 < O (3/R\). Parametrically, this is a factor A~ 2 I 3 smaller than the range of A- values in which a 
nucleon 2GPD has support. If one neglects the A-scale as being small, then the single GPDs become standard parton 
distributions and the 2 GPDs become two-parton distribution functions. In particular, the single nucleon contribution 
to the 2GPD in Q can be approximated as 

G A N (x u x 2 ,A) ~ AG N (x u x 2 ,A) -> Af N {xi,x 2 ) , (15) 

where fN{x\,x 2 ) is the standard two-parton distribution function. Similarly, the double nucleon contribution to the 
2 GPD in (|9| can be approximated as 

G A N (x 1 ,x 2 ,A)^f N (x 1 )f N (x 2 ) Ft ublc (A,-A), (16) 

where /jv(x) are standard single parton distribution functions. Simplified expressions for the double hard four jet 
cross section can then be obtained by inserting equations (151, (16) into 

• Case III 

We consider first the 4 — >■ 4 contribution to the double hard 4-jet cross section in which the double nucleon 
scattering term G 2 A of the nuclear 2GPD is paired with Gp 0ublc (xi, x 2 , A) in the proton. We obtain 

da^ II " > {x' l ,x' 2 ,Xx,x 2 ) daida-2 f d 2 A doublc , , 1 X\ f / \ t 1 

- A{A-l)— r — r / — G [x 1 ,x 2 ,A) f N {x x ) f N (x 2 )F A (A) . (17) 



dtxdt 2 dt x dt 2 J (2?r) 



Taking the nucleus large enough to ignore the nucleon size, see eq. (16), one can neglect the A dependence of 
G p (xi, x 2 , A). Then one can write 

d*i"M,4,*u**) _ A{A _ l) do^da 2 f ^A G dou W e ( ^ )4)0)/jv(xi)/iv(:C2)jFl(A) . (lg) 



dt x dt 2 " dt 1 dt 2 J (2?r) 2 

Here Gp° uhlc (x' 1 , x' 2 ,0) = / p (x' 1; x' 2 ) is the double parton distribution function, and Jm denotes standard nucleon 
pdfs. For a simple parametric estimate, the form factor F\{A) can be viewed as a step function with support 
for A 2 < 3/R\ ~ t4~ 2 / 3 . Therefore, the contribution (18) is 0(A 4 / 3 ) which makes it A 1 / 3 -enhanced compared 
to all contributions discussed in section [ill This can also be seen after Fourier transform to 6-space, if one recalls 
that T(b) cx A 1 / 3 for typical b < R A - A 173 , 

cri III \x' 1 ,x' 2 ,x 1 ,X2) = f P (x' 1 ,x' 2 ) da^jx^x^ da p 2 f ot (x' 2 ,x 2 ) (A - 1) 
dhdh dh di 2 A 




(19) 



Here we expressed the four jet cross section in term of full dijet differential cross sections, defined through hard 
parton cross sections as: 

^( x [, Xl ) = f N ( Xl )f p (x' 1 )^(x' 1 ,x 1 ). (20) 
at at 

Except for the (1 — 1/A) correction term, this form of the double hard 4-jet cross section was given first in |28j . 
We note that in evaluating J T 2 (b)d 2 b in (19), short-range nucleon interactions can be taken into account. For 
A ~ 200, the resulting corrections are on the level of a few percent |55] . 

• Case VI 

The double nucleon scattering term G 2 A N enters also in the 3 — > 4 contribution to the double hard 4-jet cross 
section, where one parton of the proton splits into two partons with momentum fractions x[, x' 2l 

^ = A(A-1)^^ /-^Gr^K,4,0)/ JV (x 1 )/iv(x 2 )Fi(A)ocA 4 / 3 . (21) 
dtidt 2 dti dt 2 J (2ir) z 
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This term has the same parametric A 4 / 3 -enhancement as (18). However, as we explain now, its relative weight 



compared to a^ 111 ^ is significantly smaller in pA than the corresponding relative weight in pp collisions. To see 
this, let us recall first that in pp collisions the 4—^4 contribution involves a A-integral over the fourth power of 
.Fjv(A). This follows for instance from writing in the mean field approximation each of the two 2 GPDs in |2]) as 
the product of two single GPDs iG(x, Q 2 , A) = /n(x, Q 2 ) F/v(A). On the other hand, in proton-proton collisions 
the A-integral of the 3 — > 4 contribution involves only two powers of F/v(A), since G p mgle (xi, x 2 , A) corresponds 
to a point-like perturbative splitting and its A-dependence is thus negligible lj. In general, since F N (A) peaks 
at A = and falls off steeply with increasing A, the A-integral over Fj^(A) is larger than that over F^(A). 
This results in a geometrical enhancement of the 3 — > 4 contribution in pp relative to the 4 — ¥ 4 contribution. 
For a numerical estimate of this enhancement in pp, one may take recourse e.g. to the F(A) — ^ | ^i^^y [36|. 
for which 



d 2 A 
J2^) 2 



d 2 A 
(2^ 



127T 2 / 28tt 2 



7/3. 



(22) 



We note that this ratio is rather robust against changes of the functional shape of F{A); for instance, an 
exponential form of F (A) would yield a ratio 2 rather than 7/3. In contrast to this geometrical enhancement 
in pp, we have found here that in pA collisions in the limit of very large A, the A-integrals in ( 18 1 and ( 21 ) are 



the same and a geometrical enhancement as in (22 1 is missing. In addition, the single parton in a 3 — > 4 process 



could belong to each of the colliding protons in pp, whereas this combinatorial factor 2 is obviously absent in 
pA. In summary, there is a geometrical enhancement factor of 3 — > 4 relative to 4 — > 4 processes in pp that 
equals 7/3x2 = 14/3 ~5 and that is clearly absent in pA collisions for sufficiently large A, 



AVI) 
Za 

Am) 



T (VI) 

Za 

Am) 



^ = const(A)| A>>1 - i 



(23) 



IV. TOTAL CROSS SECTION IN INDEPENDENT NUCLEON APPROXIMATION. 



In summary, working in the independent nucleon approximation and neglecting the nucleon size compared to the 
nuclear radius Ra, one can write the double hard four-jet cross section in pA as the sum of terms that are linea r in 
A (see discussion of a\ N in section Inj and cq. ^ ) and of double scattering terms that were considered in section III 

(see eqs.(19), (21) ). We define <jf^} t as the sum of all these contribution. Normalizing this sum by A times the 
corresponding cross section in a nucleon-nucleon collision, we obtain the 4-jet nuclear modification factor 



D 4?'ei / I I \ 

R pA (Xl,X 2 ,X 1 ,X 2 ) = 



dti dt 2 
S A-1 
A A 



T 2 {b)d 2 b 



Ad<72jet{x' l ,Xi) da 2 jet(x2,X 2 ) 

S dt\ 
Gp(x 1 , x 2 ) 



fp( X l)fp( X ' 2 ) 



where 



G p (x' 1 , x' 2 ) — f p (x' 1 ,x 2 ) + G p ln ° lc (x' 1 , x' 2 , 0) . 



(24) 



(25) 



We note that a closely related expression was obtained already in the analysis of Ref. [28] , where the 3 — > 4 contribution 



was not included and 1 /^4-corrections were neglected. It is useful to write ( 24 ) in the compact form 



^pA t {xi,x 2 ,x' 1 ,x' 2 ) 



1 + SW(A)K(x' 1 ,x' 2 ) 



(26) 



where the second term on the right hand side of (24) factorizes into a product of the effective transverse area S, a 
purely geometrical overlap factor 



W(A) = ^J d 2 bT 2 (b) 



(27) 




FIG. 3: Schematic representation of an interference contribution to the nucleus 2GPD in which the two patrons of momentum 
fractions xi, £2 are drawn from two different nucleons in the amplitude, but in which the assignment between partons and 
active nucleons is swapped in the complex conjugate amplitude. As discussed in section |Vj these contributions (cases VII and 
VIII) are expected to suppressed due to constraints on the longitudinal light-cone fractions of the two active nucleons. 



and the normalized longitudinal parton correlation function 



K(r' t' ) - gpOgij^k^ (281 



The second term SW(A)K(x' 1 ,x' 2 ) on the right hand side of (26) corresponds to the parametrically j4 4 / 3 -enhanced 
contributions to the 4-jet double scattering cross section that we have discussed in section |III| For a lead nucleus 
with Wood-Saxon density profile and S = 15 mb, one finds SW(A) ~ 2.1 [37]. This shows that these parametrically 
enhanced terms give indeed the largest contribution for sufficiently heavy nuclei. However, terms with linear A- 
dependence must not be neglected for numerical estimates, since they constitute about 1/3 of the four-jet cross 
section. 



Equation (24) or (26) summarizes one of the main results of this paper. It demonstrates that any x\- and x' 2 



dependence of the nuclear modification factor provides information about the longitudinal correlation ( 25 ) of two 
partons in the nucleon, since it constrains directly the ratio K. It is in this sense that the nucleus provides a non- 
trivial filter for analyzing the multi-parton structure of the proton. We further note that the effective area S is 
operationally defined as the ratio of the double hard 4-jet cross section and the product of two dijet cross sections in 
pp, and thus it can be a function of S = S(x±,X2, x[, x 2 ). However, 

K( Xl ,x 2 )- SW{A) (29) 



can depend only on the momentum fractions x' x , x' 2 in the proton. Any deviation of the ratio ( 29 ) from unity would 
be an unambiguous signal of longitudinal momentum correlations in the proton. We note that the numerical analysis 
of [3J, based on dominance of 3 — > 4 processes, suggests K ~ 1.2 in the kinematic region under consideration. Vice 
versa, since the 3 — > 4 contributions involve partonic 1 — > 2 splittings in the nucleon and thus lead dynamically 
to longitudinal momentum correlations, any tight bound on \K(x' 1 ,x' 2 ) — 1| will put significant constraints on the 
role of 3 — >• 4 processes. In the extreme case when K(x' 1 ,x' 2 ) = 1, one would have to conclude that the dynamical 
origin of the anomalously small effective area 5 is purely transverse. We argue that these considerations motivate 
an experimental study of the nuclear modification factor (|24|) and the corresponding longitudinal correlation function 



(29) in the upcoming pA run at the LHC. 

We note as an aside that the corresponding result for proton-deuteron scattering can be calculated, based on the 
explicit form of the two-body deuteron form factor F^ ouble (A, — A). Since the nucleons in the deuteron are strongly 
correlated in the center of mass frame, i^ oublc is simply expressed through the one body deuteron form factor Fd as 
F5 (A, —A) = Fd(4A 2 ), similar to the case of the Glauber scattering, cf. [Mj. Inserting this expression into (21 ) 



and performing the A-integral, one obtains ~ 1 + I.Q7K for x t = 0(0.01). We thank V. Guzey for performing 

this numerical integration and communicating the result. The case of pD scattering was considered recently in |30j 
in coordinate space formalism, though no numerical results were reported. 
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V. DOUBLE SCATTERING - INTERFERENCE TERM (CASES VII AND VIII) 

Our discussion so far was based on a picture in which the nucleus is viewed as an independent superposition of 
nucleons, even if the parton distribution functions inside each nucleon may differ from those in a free nucleon. If 
one does not assume this picture, further interference contributions to the double hard four-jet cross section are 
conceivable. Here, we discuss the form of these contributions for completeness. We note at the beginning that we do 
not have a complete framework for calculating them, but we shall give some arguments of why we expect them to 



provide at best very small corrections to the expression (24). 



As depicted in Fig. 3, it is possible to write down a diagram where in the initial state a parton with x\ from the 
nucleon " 1" is involved in the two to two process while this nucleon absorbs the parton with x 2 m the final state. 

il f 

This process changes the longitudinal light cone fractions a^ 2 °^ ^ ne ^ w0 ac tive nucleons between the initial and final 
state, 

a\ — x\ + x 2 — a{; a 2 + x\ — x 2 — ot 2 . (30) 

It thus implies that the momentum transfer A now has also a nonzero longitudinal component, which in the nonrel- 
ativistic approximation can be written as 

A = (A 3 = (xi - x 2 )m Nl A t ). (31) 
Consequently, in close analogy to the discussion of ( fl4| ), the two-nucleon form factor takes now the form 

F A {A, —A) » Fl((xi - x 2 ) 2 m 2 N + A 2 ) a exp(-((a; 1 - x 2 fm 2 N + A 2 )R 2 A /S) (32) 



in the mean field approximation. Here, the transverse factor is of the form of (14 1. The additional suppression factor 
exp(— {x\ — x 2 ) 2 m 2 N ■ R\/3) arises for significant differences in the longitudinal momentum fractions. For a typical 
value of the nuclear radius at large A, R A ~ 6 fm, one finds a strong suppression in the range \x\ — x 2 \ > 0.03. 
For \x\ — x 2 \ < 0.03, this suppression factor is less important and it vanishes for x% — x 2 . Naively, one may expect 
that in this phase space region such a contribution is thus enhanced. On general grounds, however, one may argue 
that this scattering diagram is also strongly suppressed in the range \xi — x 2 \ < 0.03, since a momentum exchange 
between nucleons should imply that the active nucleons in the nucleus are close in configuration space. Here, we do 



not attempt to estimate such potential further suppression effects further. We merely observe from (32) that such an 



additional contribution, if not negligible, will correct our main result (26) 



R 4 v >f{x 1 ,x 2 ,x' 1 ,x' 2 ) = l + SW(A)K{x , 1 ,x' 2 ) (1 + F corr ) , (33) 



where the correction factor F corr is expected to be small. In the case considered in (|32|), this correction factor will 
take the explicit form i^orr = const • exp(— {x\ — x 2 ) 2 m 2 N R 2 A /i), where the proporti 
unity. 



ity factor is much smaller than 



VI. TWO NUCLEON - TO - ONE NUCLEON INTERFERENCE (CASE IX) 

We finally consider the process where the two partons active in the double hard scattering are drawn from the 
same nucleon (with light-cone momentum fraction a\) in the amplitude, but from two different nucleons (with light 
cone momentum fractions a^, a^) in the complex conjugate amplitude, see Fig|4j The momentum fraction of the 
spectator nucleon in the amplitude is denoted by a 2 . Longitudinal momentum conservation implies constraints on 
such processes. In particular, a 2 + {3 = and x 2 < ft. Since the light cone momentum fractions on deviate from unity 
only by Fermi motion, this will constrain contributions of the diagram Figj4] to relatively small momentum fractions 
x 2 . In addition, we expect that exchanges as depicted in Figj4] can arise only between nucleons neighboring in impact 
parameter so that this diagram should be also suppressed by geometrical factors. Qualitatively, the limitations to 
sufficiently small x and a further suppression of geometrical origin are the same two suppression factors that we 
identified for the cases VII and VIII discussed above. To become more quantitative is difficult and lies outside the 
scope of the present work. We note, however, that diagram Figj4] is closely related to the leading twist shadowing 



correction, discussed in Ref. [35]. We now turn to a more general discussion of the corrections to (24) that may arise 
from the nuclear dependence of parton distribution functions. 
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FIG. 4: Schematic representation of an interference contribution to the nucleus 2GPD in which the two partons are drawn from 
a single nucleon in the amplitude while the parton '2' is associated to a different nucleon in the complex conjugate amplitude. 
This process requires the exchange of additional constituents between both active nucleons, as indicated by the grey blob. 



VII. NUCLEAR DEPENDENCE OF SINGLE PARTON DISTRIBUTION FUNCTIONS 



In our discussion so far, we have neglected nuclear modifications of the single parton distribution functions. Standard 
parametrizations of these modifications are based on a linear relation between the parton distributions fi/ p (x, Q 2 ) in 
a proton and the parton distributions fi/A(x, Q 2 ) per nucleon in a nucleus [3i 



h/Afa Q 2 ) = Rf& Q 2 ) fi/ P {x, Q 2 ) . 



(34) 

In all parametrizations that are consistent with linear Q 2 -evolution [3"B1 BOTES] , the npdf factors Rf(x,Q 2 ) at fixed 
x approach unity with increasing Q 2 . Since the jet production considered here is a hard process with typically 
Q 2 » 100 GeV 2 , we expect \Rf (x,Q 2 ) - 1| < 3% over the entire x, and Q -range relevant for four-jet production at 
the LHC. This assumption can be verified experimentally by checking that the nuclear modification factor for dijet 
production, 



p2jct 
n pA 



d<J lfet 

dt\ ' 



A 



da p Z 



2 jet 



du 



(35) 



deviates from unity by less than 3% in the kinematical range used to measure the four-jet cross sections. An 0(3%) 
uncertainty from "npdf corrections" to the dijet cross section entering the norm in (24) is expected to result in an 
0(6%) uncertainty on the nuclear modification factor for double hard four-jet production. The central question is 



whether a correction to ( 24 ) of this order is small enough to be neglected in the analysis of the longitudinal two-parton 
correlation function K{x\, x 2 ). This depends on how the size of the deviation of K(x\,x 2 ) from unity compares to 
the size of the npdf correction of (24 1, and we shall distinguish below two different cases. 

Let us discuss, however, first how the nuclear modification of longitudinal parton momenta in the nucleus can 
be taken into account in the formulation of 2GPDS. To this end, one needs to specify correction factors for the 
three contributions to the nuclear 2GPD listed on the right hand side of equation Q. In principle, this requires more 
information than what is contained in the npdf- fits based on ( 34 ) . Adopting the picture of the nucleus as an incoherent 



superposition of nucleons that have parton distributions shifted according to (34), one arrives at corrections of the 
simple form 



G A N (x 1 ,x 2 ,0) 

G double, IN / n \ 
A (xi,x 2 ,0) 



R?(x 1: Q 2 ) Rf(x 2 ,Q 2 ) G 2 A N ( Xl ,x 2 , 0) 

-rdouble,lAf 



Rf(x 1 ,Q 2 )Rf(x 2 ,Q')G° J 



(xi,x 2 ,0) . 



(36) 
(37) 



In principle, nuclear pdfs depend also on the impact parameter b (see Rcf. [38, 39 for first parametrizations), and 
this 6-dependence could be accommodated in (36) and (37). This, however, will be a small correction on top of the 



correction discussed here, and we neglect it in the following. 
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The third contribution Q 8 ^ 3 e ' 1N ^ the nuclear 2GPD does not have an npdf correction factor of this form. Rather, 
since the arguments X\, x 2 in Q 8 ™s le ' 1N result dynamically from the splitting of a single parton, the npdf corrections 
will be determined by an integration over the available phase space that may be written formally as 



GT 9le ' 1N (x 1 ,x 2) 0)^O (Rf(x 1 + x 2 ,Q 2 )) GT 9le ' 1N (x 1>X2 , 0) . 



(38) 



To determine npdf corrections to ( |24| in the most general case, one would have to specify (38) fully and then repeat 
the calculations in sections ITflTV based on equations p6|)-(38|. Because of the more complicated form of (38), the 



explicit results for the npdf corrections to (24) would be relatively involved. Here, we restrict our discussion by 
considering two limiting cases: 

First, we observe that if (38) makes a numerically important contribution to the total nuclear 2GPD Q, then 
longitudinal two-parton correlations in the nucleus are expected to be large. This is so, since (38) is dynamically 



generated by a perturbative parton branching that inevitably results in significant longitudinal correlations. Indeed, 
according to the model of [3] one can expect a ~ 20% deviation of K(xi,x 2 ) from unity if the double hard four jet cross 
section is dominated by 3 — > 4 processes. Therefore, if the term (38) is relevant, we expect that K(x\,x 2 ) — 1 > 0.1 
and that the npdf- correct ions to (24) are small. In this case, K[xi,x 2 ) — 1 is much larger than the expected npdf 



corrections and it can thus be extracted safely from (29) without taking npdf corrections into account. We emphasize 



that the validity of this procedure can be checked experimentally by measuring the nuclear modification factor for 
dijet production in pA. 

Alternatively, if an experimental determination of K(xi,x 2 ) — 1 from (29) yields values K(xi,x 2 ) — 1 <C 0.1, the 
contribution of the 1 — > 2 splitting term (38) to the nuclear 2GPD can be expected to be small, and one can justify 



hence the use of a simplified npdf correction of the form 2 Ga(xi, x 2 ) — > Rf(xi,Q 2 )Rf(x 2 ,Q 2 ) 2 GA{xi,x 2 ). In this 



case, the npdf corrections to the contributions ([7]), (JT9j) and (21) amount to multiplying all three expressions with 
two powers of the nuclear modification factor for dijets, R 2 J A e . As a consequence, the right hand side of (24 1 gets 



multiplied by two powers of RJa , and the npdf correction to ( 29 ) takes the form 



corrected 



x 2 ,xi,4)/(i^ 2 



SW(A) 



(39) 



Experimentally, comparing the values obtained from ( 29 ) and ( 39 ) provides a direct way of estimating the importance 



of npdf corrections on the interpretation of if as a normalized longitudinal two-parton correlation function. 



VIII. CONCLUSIONS. 



In summary, by analyzing 2GPDS for the nucleus in the many nucleon approximation, we have derived a compact 



expression (24) for the nuclear modification factor of the double-hard four jet cross section in pA collisions. Based on 



this main result, we have outlined an experimental strategy for determining the normalized longitudinal two-parton 
correlation function K(xi, x 2 ) in the proton by combining data from pp and pA collisions. Finally, we have discussed 
how nuclear modifications of parton distribution functions can be taken into account in this analysis. 

In early 2013, the LHC is scheduled for a 4-week-long proton-nucleus run. The main motivation for this pA program 
at the LHC is to constrain the parton distributions in the nucleus and to provide important benchmark measurements 
for the LHC heavy ion programme. As illustrated by the calculation of K{x'^x'^) in this paper, however, proton- 
nucleus collisions may also contribute to further constrain the multi-parton structure of the proton, thus probing the 
proton in nuclear collisions rather than probing the nucleus in collisions with protons. 
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